Lower bounds on entanglement of formation for general Gaussian states 
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We derive two lower bounds on entanglement of formation for arbitrary mixed Gaussian states 
by two distinct methods. To achieve the first one we use a local measurement procedure that 
symmetrizes a general Gaussian state and the fact that entanglement cannot increase under local 
operations and classical communications. The second one is obtained via a generalization to mixed 
states of an interesting result already known for pure states, which says that squeezed states are 
those that, for a fixed amount of entanglement, maximize Einstein-Podolsky- Rosen-like correlations. 
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I. INTRODUCTION 

The quantification of the amount of entanglement a 
quantum system possesses is still an open problem in 
Quantum Information Theory. Restricting our attention 
to bipartite systems, i. e., systems composed of two sub- 
systems, we have one measure of entanglement, entangle- 
ment of formation (EoF) 1], which has a clear physical 
meaning. Given an entangled state p, the EoF for this 
state expresses the number of maximally entangled states 
we need to create p . The formal definition of the EoF 
is: 



E f (p) = Mj2PjEty j ), 



(1) 



where we take the infimum over all pure-state decompo- 
sitions of p = J2jPj \ipj) (i>j\> J2jPj = 1 and E(ipj) is 
the von Neumann entropy of the pure state ipj . 

The analytical minimization of Eq. is not an easy 
task. Dealing with two-qubit systems, which are the sim- 
plest entangled bipartite systems, Wootters obtained 
an analytical expression for the EoF and Giedke et al Q| 
derived an analytical expression for the EoF for symmet- 
ric Gaussian states. 

Gaussian states are very useful in quantum-optical im- 
plementation of several quantum information protocols. 
(Quantum cryptography 5] is an important example.) 
Hence, a complete characterization of the amount of en- 
tanglement of Gaussian states is desirable. The natural 
next step is the search for an analytical expression for 
the EoF for arbitrary Gaussian states. 

In this article we give two analytical expressions that 
furnish lower bounds for the EoF for Gaussian states. 
We employ two different methods to derive such lower 
bounds. The first lower bound is obtained using a lo- 
cal measurement procedure derived by Giedke et al [|J 



which symmetrizes a general Gaussian state and the fact 
that entanglement cannot increase under local operations 
and classical communications (LOCC). The second one 
is derived via a generalization to mixed states of an inter- 
esting result derived by Giedke et al |4|, who show that 
squeezed states are those that, for a fixed amount of en- 
tanglement, maximize Einstein-Podolsky- Rosen- like cor- 
relations. These lower bounds are also useful to rule out 
several possible candidates for the analytical expression 
of the EoF for arbitrary Gaussian states, as we illustrate 
in this article. 



II. FIRST LOWER BOUND 

Let us begin setting the notation used in this article 
and some properties of Gaussian states. Consider a bi- 
partite Gaussian system p of two modes described by 
the annihilation operators a,j — (Xj + iPj)/y/2, where 
j = 1,2 and [Xj,Pj>] = iSjj'- This system can be alter- 
natively described by its characteristic function 

X (r)=tx[pD(r)], (2) 

where r — (xi,p\,X2,P2) T is a column real vector and 

D(r) = e ~ i ( XlXl +P lPl + X2X2 +P 2P2 ) (3) 

Eq. (J2J uniquely defines a state p and for Gaussian states 
it can always be put in the following form: 



x (r) = e-V ^~ ld r , 



(4) 
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where T means transposition, 7 is a 4 x 4 matrix which is 
called correlation matrix (CM) and d is a 4— dimensional 
real vector. The first moments of a Gaussian state (Xj) 
and (Pj) can always be set to zero using local unitary op- 
erations, which implies that we can work with zero mean 
Gaussian states when studying entanglement properties 
of such systems. The matrix elements 7^ of the CM can 
be calculated directly from the density matrix p by the 
following formula: 

7„- = tr [(RiRj + R 3 Ri)p] - 2tr[i? l p]ir[i? J p], (5) 



2 



where R = (X\, Pi, X2, Pz) T ■ A matrix 7 represents a 
realizable physical state iff it is strictly positive, real, 
symmetric and satisfies 0: 



7 > j t i~ x j, 



(0) 



where J = ©t =1 Ji is a 4 x 4 matrix with Ji 

-1 

1 

A Gaussian system can also be represented by its 
Wigner distribution W(r). Assuming that we are work- 
ing with a zero mean Gaussian state we have Q: 



W(r) = 



1 1 



-r 'yw r 



7T 2 y/detjw 



(7) 



The CM's 7 and jw are related by the following relation: 

(8) 



lw = J T 1 1 J- 



These two CM's can be brought to the following stan- 
dard form by suitable local symplectic transformations 



6] 



A C 
C T B I ' 



(9) 



where 
A = 






, B = 


' m 


° 









n J 




m J 



c = 



k x 
k v 



The same set of equations apply to jw- 



Aw Cw 
<^w B w 



where 



, N \ D (MO 
Aw ( N I , B w = I Q M 



w 



K x 

K„ 



(10) 



(11) 



(12) 



(13) 



The four real parameters (n, m, k x , k p ) completely char- 
acterize a two mode Gaussian state and they are related 
to the four local symplectic transformation invariants as 
follows 0: 



I\ = n = vcletA, 



Jo = m = vdetB, 



(14a) 

(14b) 
(14c) 



h = k x k p = detC, 
h = nm{k 2 x + k 2 p ) = tr {AJ^ C j{ BJ^ C T Ji) . (14d) 

Alternatively the four real parameters (N,M,K x ,K p ) 
also completely specify a two mode Gaussian system. 



They can be also obtained by local symplectic trans- 
formation invariants. These invariants, which we call 
Wi,W2,Ws and W4, satisfy Eq. (fPp . where we change 
A, B and C by Aw,Bw and Cw and {n,m,k x ,k p ) by 
(N,M,K x ,K p ). 

We now pass to the derivation of the first lower bound. 
A symmetric Gaussian entangled state a is completely 
specified by its CM (see Eq. where n = rn = h. 

(From now on, every parameter associated with a sym- 
metric Gaussian state will be represented by a tilde on 
top of it.) Let us assume, without loss of generality, 
k x > and k p < . The EoF for this symmetric state 
is H: 



Ef(a)=f 



(h - \k x \)(n - \k p \) 



where, 

f(S) = c+(S) log a M*)] - c-(S) log a [c_(*)]. 



(15) 



(16) 



Herec± = {5-^ 2 ±S^ 2 ) 2 /A. Using Eq. we can write 
the EoF given by Eq. I)15|l in terms of invariants: 



Ef(a) = f 



2hh 



(17) 



Using Eqs. JSJ) and l|14|l we obtain the following relations 
among the invariants of the 7 and jw matrices: 



W2 

w 5 



Wi 

W 5 



W3 
w 5 



1 

w 5 



(18) 

where W5 — detjw and I5 — detj. 

Therefore, due to Eq. I|18|l the EoF for our symmetric 
Gaussian state, Eq. (|17fl . can be expressed as: 



Ef(a) - / 



'Wi -W 3 - VW 4 - 2WxW 3 



w 5 



(19) 



But Giedke et al |6( have shown that a general bipartite 
Gaussian system p can be transformed to a symmetrical 
bipartite Gaussian system a using LOCC. This implies 
that Ef(p) > Ef(a). Schematically we have: 



LOCC 

p — > a 



E f (p)>E f (a). 



(20) 



Our only task now is to rewrite Eq. (|19|) in terms of the 
invariants of the 7 matrix of p. It is in this step that we 
use Giedke's symmetrization procedure. 

Given a general bipartite Gaussian system p and its 
jw matrix, where we assume, without loss of generality 
that N > M, we can achieve by means of local operations 
a symmetric state with the following 'jw matrix : 



lw 



A w 



c 



w 



C w Bw 



LOCC _ 
, — > lw = 



A w 

IT 



C 



w 



C w B w 



, (21) 
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where 
A w = 



N < 



cos 2 6+M sin 2 ! 





jVcos 2 6+NM sin 2 6 
cos 2 6+M sin 2 



(22) 



e = W 5 + VWi( V^VW^ - #2) tan 2 9, 



(28) 



D TIr _ cos 2 6+M sin 2 6» " I /OQ^i 

^ _ 1 sin^ + Mcos 2 ^ 1 ' l " ! 



a 



w 



Kg cos g 
s 2 0+Msin 2 e 



iCp cos 9 



(24) 



/3 = IU 4 — 2W2W3 + tan 2 6> (W 4 - 2W 3 - W|) VW2 



+ (1 - iy 2 )X 2 v/Wi 



(29) 



tan" 4 6< = 



7V 2 -M 2 



M - 7V(7VM - ■ 



(25) 



Eq. (|2*5|l guarantees that detA^ = detB^v. This condi- 
tion is the statement that the Gaussian system with the 
Tiv above is symmetrical 0. 

Using Eqs. fTJIZ232i and the assumption that \K X \ > 
\K P \ Q we can write Eq. (|19J) as follows: 



where 



Ef(cr) = I 



W 2 -W 3 



-VP 



'W 2 tan 2 9, 



(26) 



(27) 



tan^ 9 = 



Wx - W 2 



Kt 



- VW{^fWW2 - Kl) ' 



w 4 + y/wl - m x w 2 wl 
2VW1W2 



(30) 



(31) 



Now Using Eq. (|18|l we can put Eq. (|26[) in terms of the 
invariants of the 7 matrix. Hence, if we work with 7 in 
its standard form given by Eq. @, where we assume, 
without loss of generality, that \k x \ > \k p \, Eq. ffify is 
rewritten after a tedious but straightforward algebraic 
manipulation as [To| : 



Ef(a)=f 



I nmh(n, m) — k x k p h(m, n) + \mk x — nk p \Vh(n, m)h(m, n) 

g(n, m) 



(32) 



where 

h(n, m) = n — m(nm — k p ) (33a) 
g[n,m) = m(l — m 2 ) + nfc 2 . (33b) 

Eq. llol'l) is our first lower bound for the EoF for gen- 
eral Gaussian states. It is worthy noting that this lower 
bound reduces to the EoF for symmetric Gaussian states 
whenever n = m. 



III. SECOND LOWER BOUND 

A two-mode squeezed state [4] is a symmetric Gaussian 
pure state that belongs to the Hilbert space H = Tti®H 2 



and is described by the following vector: 

= vTT E tanh » l">i ® l n >2 > ( 34 ) 

cosh r 

v ; n=0 

where |n)^ is the n-th Fock state, that is, a^aj \n) ■ 
n \n)j, j=l,2 and r E (0, 00) is the squeezing parameter. 

There exists an interesting relation between squeezed 
states and EPR-correlations, which Giedke et al 
proved in their proposition 1: Given a squeezed state 
|* s (r)) and an arbitrary pure two-mode state \ip) then, 
if they have the same EPR-correlation, the squeezed state 
is the least entangled. In other words, if we call A($ s (r)) 
and A(-0) the EPR-correlations for the two mentioned 
states and if A(* s (r)) = A(ip) then, E(ip) > E(^ s (r)). 
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The EPR-correlation is defined as £| : 



A(^)=min<jl,-[A|(X 1 



X 2 



A 2 (P 1 + P 2 )] 



(35) 

where A^(Rj) = (Rj)^ — {Rj)^ is the dispersion of the 
observable Rj . The above expression measures the degree 
of non-local correlations, and is zero for the original EPR- 
state (3, El ■ This means that the more a system is non- 
local the more Eq. approaches zero. We say that a 
system with the minimal A(ip) has the maximal EPR- 
correlation. For our squeezed state the EPR-correlation 
is H: 



A[* a (r) 



-2r 



(36) 



The EoF, which is equal to the von Neumann entropy 
for the squeezed state is [4(: 

E[^ s (r)} = cosh 2 (r)log 2 [cosh 2 (r)] 

— sinh 2 (r) log 2 [sinh 2 (r)] . 

And it is shown that Q : 

E[* a (r)] = /(A[*.(r)]). 

Giedke et al Q have shown that / : (0, 1] — * [0, oo) is a 
convex and decreasing function of its argument 
as Eq. (|36[) can have any value between zero and one, the 
EoF for a squeezed state can assume any value between 
zero and infinity. This property of the EoF for squeezed 
states, i. e., that they can assume any value, is an essen- 
tial ingredient in our generalization of Giedke's et al |J| 
proposition 1. Let us now state and then prove the fol- 
lowing theorem which is a generalization to mixed states 
of Giedke's et al Q proposition 1. 

Theorem 1 For all bipartite Gaussian systems p we 
have Ef(p) > Ef(a), if A(p) — A(a) and a is a sym- 
metric Gaussian mixed state. 

Here A(p) is analogously defined as in Eq. (|35|l . 
Proof: Applying a suitable symplectic local transforma- 
tion in the standard form of the 7 matrix of a |4j, [12J we 
see that the EPR-correlation for this transformed matrix 

But the amount of 



is A(<t) = J [n - \k x \j [n - \k 

entanglement is invariant by local symplectic transfor- 
mations. This means that Ef(a) = f [A(er)] = / [A(p)]. 
Let us write p as 



P = J2Pj \<Pj) {fj 



(39) 



where the above decomposition is the one that furnishes 
the EoF of p, i. e., 



(40) 



Using the above expansion of p we have that 



Ef{ci) = f 



<f 



<X>/[A(¥>j 



(41) 



The first inequality is a consequence of the concavity of 
A(p) (see Appendix) and the fact that / is a decreasing 
function of its argument Mj. The second inequality is 
due to the convexity of / Q . We now use the fact that 
a squeezed state can assume any value of entanglement. 
For each pure state in the decomposition of p above we 
associate a squeezed state with the same amount of en- 
tanglement: E(ifj) — E[$> s (rj)]. Therefore we have the 
following relation for the EoF of p: 



= 5>/[A[* s (r,)]]. (42) 

j 

Now due to the proposition 1 of Giedke et al 3 we 
(37) know that A((fj) > A[H> s (rj)]. Hence, using this fact 
in Eq. lEl) and that / is a decreasing function of its 
argument we have: 



(38) 



Hence, Ef{o) < ^Tpjf [A(^-)] < ^^Pjf [A [* a (rj)]] . (43) 



Combining Eqs. (|42|l and (|43|) we see that 



E f (a) < E f (p). D 



(44) 



The above theorem tells us that for mixed states the 
symmetric states are those with less EoF given an EPR- 
correlation. It is interesting to note that a can be any 
symmetric state, including symmetric states written as 
superpositions of squeezed states. 

The previous theorem automatically gives us a lower 
bound for the EoF for general Gaussian states. Using 
Eq. ISI) we get: 



E f (p)>E f (a) = f[A(p)}. 



(45) 



We now implement a local symplectic transformation in 
the 7 matrix of p, Eq. ©, before calculating the EPR- 
correlation. (It does not alter the amount of entangle- 
ment, since it is equivalent to a unitary local transforma- 
tion in the density matrix p.) This transformation can 
be viewed as an extension to non-symmetrical Gaussian 
states of the transformation introduced by Giedke et al 
for symmetric states. This transformation multiplies 
Xj by [(n + m)/2 - |fc p |]/[(n + m)/2 - |fc x |]V4. P 3 is 
divided by the same quantity. Now calculating A(p) we 
get the following expression for our second lower bound: 
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Ef(p)>f 




(46) 



Again we see that this lower bound reduces to the EoF 
for symmetric systems whenever n — m. It is important 

to note that when 



n+m. 
2 



\K\) - > 1 we 

have A(p) = 1. For such cases this lower bound is not 
useful since it simply shows that Ef(p) > 0. 



state. Consider, just for illustration, the functions 



h = I 



(48) 



IV. DISCUSSION AND EXAMPLES 

We now employ the two lower bounds derived previ- 
ously, Eqs. i|32|) and (|46[1 . to see their usefulness in analyz- 
ing some Gaussian states. For completeness we present 
in terms of the invariants (n,m,k x ,k p ) three inequali- 
ties they must satisfy to be considered parameters that 
describe physically realizable entangled Gaussian states 
0. We will assume, without loss of generality m > n 
and \k x \ > \k p \. 



det7 + 1 > n 2 
nm — k 2 > 1 
det7 + 1 < n 2 



(47a) 
(47b) 
(47c) 



The last inequality is the restriction a 7 matrix must 
satisfy to represent an entangled Gaussian system. 

The table below shows six entangled Gaussian sys- 
tems and the values of their two lower bounds (LBl and 
LB2). These six Gaussian systems are very representa- 



TABLE I: The first column shows the parameters of the 7 
matrix when written in its standard form. The second and 
third column represent the two lower bounds for the EoF for 
mixed Gaussian states. Lower bound 1 is given by Eq. (I32H 
and lower bound 2 is given by Eq. 14611 . 



Thf TTi; kx 5 kp 

1.5, 2, 1.2, -1 
1.5, 2, 1, -1 
2, 3, 1.8, -1.2 
1.7, 2.6, 1.3, -0.9 
2, 3, 1.7, -1.2 
2, 2.5, 1.3, -1.2 



LBl 
0.14635 
0.08687 
0.02448 
0.00549 
0.00725 
0.00173 



LB2 
0.28919 
0.14672 
0.00681 


0.00142 
0.00001 



h = I 



\ 



-\k 




\k T , 



(49) 

Both fi and f 2 reduce to the EoF for symmetric states 
when n = m. For the Gaussian states with (n, m, k xi k p ) 
= (2,2.5,1.3,-1.2) we have LBl = 0.00173 > /1 = 
0.00091 and for (n,m,k x ,k p ) = (1.5,2,1.1,-1) we get 
LB2 = 0.208853 > f 2 = 0.18621. These results show 
that fi and /a cannot be proved to be the EoF for gen- 
eral Gaussian systems since we have lower bounds for the 
EoF that are greater than fi and fi- 



V. CONCLUSION 

We presented in this article two lower bounds for the 
EoF of general Gaussian two-mode systems. They were 
obtained by two distinct methods. 

The first lower bound, Eq. (32), was derived using an 
interesting procedure derived by Giedke et al [6j that 
symmetrizes by local operations a Gaussian state and 
the well known fact that entanglement does not increase 
under LOCC. 

The second lower bound, Eq. (|46() . is a corollary of 
theorem 1, which can be interpreted as an extension of 
a previous result obtained by Giedke et al : given two 
pure bipartite systems with the same amount of entan- 
glement, the squeezed states are those with the maximal 
EPR-corrclation. Our theorem generalizes this fact to 
mixed states in the sense that symmetric Gaussian states 
are shown to be states with maximal EPR-correlation for 
a fixed amount of entanglement. 



tive. Looking at their lower bounds we see that depend- 
ing on the parameters of the system LBl or LB2 is the 
strongest lower bound. For example, the first two Gaus- 
sian systems have LB2 as the strongest lower bound but 
the four last Gaussian systems have LBl as the strongest 
one. LBl and LB2 are also useful for discarding possi- 
ble candidates for the EoF of a general mixed Gaussian 
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APPENDIX: PROOF OF CONCAVITY OF A(p) obtain for the given expansion of p and for J2j Pj^i^j) 

the following expressions: 

We need to prove that A(p) > J2j Pj^(4>j), where 
p = ^2jPj \4>j) {4>j\- Applying the definition of A(p) we 



A(p) = min { 1, - 



E« {( x % + ( p %J - E« - E* <^ 



(A.l) 



3 j 



(A.2) 
(A.3) 



where X = X\ — X2 and P = Pi + P%. The inequality 
is a consequence of the fact that we may have at least 

one ( X2 ) h + ( p2 )^ - ( X )% - (P)% > 2 - Lookin S at 
Eq. 1A.2|1 we see that it is not greater than 1. Thus, if 



Eq. I)A.1(I is equal to 1 we see that A(p) > J2j Pj^i'Pj)- 
But if it is less than 1, A(p) > J2j Pj^i'Pj) if the follow- 
ing inequality is satisfied: 



E^wJ + (e <52 Pj [(x)l + (p); 



(A.4) 



Applying the Cauchy-Schwarz inequality 0] for an ob- Eq. i|A.4|) is always satisfied. □ 
servableflweget^j-PiW^ > (j2 3 Pj (P)^) ■ Hence, 
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